Spin-Josephson effect in antiferromagnetic tunnel junctions 
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We study Josephson-like phenomena in a tunnel junction between two itinerant antilerromagnets. 
We find at the mean-field level an equilibrium current of the staggered magnetic moment through 
the junction that is proportional to the normal state conductance and to SlxSh where Sl and 
Sr are the staggered magnetic moments on either sides. Microscopically, this effect comes from 
the coherent tunneling of spin-one charge-zero particle-hole pairs with a net wave vector equal to 
the antiferromagnetic one. We explain similarities and differences with the standard DC and AC 
Josephson effects. 



One of the most striking manifestations of supercon- 
ductivity is the Josephson effect, which consists in the co- 
herent tunneling of Cooper pairs across a junction when 
both sides are superconductors. Interestingly, this oc- 
curs even though the Hamiltonian contains only single- 
electron tunneling. The resulting current is proportional 
to sm((f), where ip is the difference between the phases of 
the two superconducting order parameters. In the pres- 
ence of a constant potential difference, gauge invariance 
requires the phase difference to increase linearly with 
time, leading to an alternating current whose frequency 
depends only on universal physical constants and not on 
any material parameter. 

Superconductivity is just one example of coherence and 
spontaneous symmetry breaking. The order parameter is 
a complex number so the broken symmetry is U(l). The 
superconducting state selects a phase and the Josephson 
effect arises from the tendency to make the phase uniform 
across the tunnel junction. 

Since there are many other types of order and corre- 
sponding broken symmetries, the question of the analog 
of the Josephson effect in such cases arises naturally [l|. 
Indeed, one should expect differences in order parame- 
ters across a junction to lead to a coherent tunneling of 
the condensed objects that exist in the broken symmetry 
state. This possibility is especially relevant in the context 
where junctions between magnetic materials are of ut- 
most importance for spintronics. In fact, theoretical pre- 
dictions have been made recently concerning the possible 
existence of a Josephson-like equilibrium spin current in 
ferromagnetic (FM) tunneliunctions in analogy with su- 
perconducting junctions 0, FM long range order 
is a realization of spontaneous S*0(3) symmetry break- 
ing and an equilibrium spin current would result from 
the exchange coupling between the magnetic moments in 
the two leads, which favors alignment of order parame- 
ters. From a Ginzburg-Landau point of view, there is, in 
the functional, a term proportional to Mi ■ where 
Ml and Mr are, respectively, the magnetic moments on 
the left and on the right of the junction. The Heisen- 
berg equations of motion thus lead to ~ Ml x Mr, 
corresponding to a spin current. 

Josephson-like phenomena between antifcrromagnets 



(AF) are interesting for several reasons. In the context 
of spintronics, it has been shown experimentally that the 
absence of a net angular momentum in AF results in 
orders of magnitude faster spin dynamics than in FM, 
which could expand the now limited set of applications 
for AF materials More generally, in AF the Neel 
order parameter breaks both lattice translation symme- 
try and 50(3) spin rotation symmetry so the situation 
is less straightforward than in FM from a mathematical 
point of view. In addition, AF are often close to super- 
conducting phases, as is found in heavy fermions, high- 
temperature superconductors and layered organic super- 
conductors. Understanding Josephson-like phenomena in 
AF is the first step towards more general studies with co- 
existing antiferromagnetic and superconducting order pa- 
rameters. Generalized Josephson effects may help iden- 
tify homogeneous coexistence in real materials. Studies 
along these lines have recently appeared for ferromag- 
netism coexisting with other types of order 

In this paper, we consider a microscopic model for sin- 
gle electron tunneling between itinerant AF. This mi- 
croscopic calculation leads to an explicit expression for 
the analog of the critical current and its temperature de- 
pendence. In addition, we show that (a) Cooper-pair 
tunneling is replaced by tunneling of a spin-one neutral 
particle-hole pair (b) time dependence introduced by ex- 
ternal magnetic fields resemble the AC Josephson effect 
but there are many differences because of the non-abelian 
nature of the problem and because spins do not couple 
to the gauge field but direcly to the magnetic field. Fi- 
nally, we briefly discuss how such effects can in principle 
be observed experimentally. 

The model The Hamiltonian for a tunneling junction 
consisting of two leads of an AF material and an insulat- 
ing barrier between them reads 

H = H L(cL,c ka ) + H R {d\ a ,d^) + H T , (1) 

where Hl(r) is the Hamiltonian of the left (right) AF, 
Ht is the tunneling part connecting the two leads, and 
c kcr( Ckcr ) ano - ^qo-^qo-) are the fermion creation (annihi- 
lation) operators of the left and right leads, respectively. 
In the following discussion the quantum numbers k and 
q will also denote implicitly the left and right lead. 
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We model the AF on each side of the junction by a 
one-band Hubbard Hamiltonian treated in the Hartree- 
Fock approximation for a static spin-density wave (SDW) 
with wave- vector Q. Without loss of generality, we as- 
sume the SDW mean field to be polarized along the spin 
quantization axis. Following Ref. [H, we write 

H L = 22 EfccLckc -—^2 c£+Qaff*SCk0 , (2) 
ka ka/3 

where efe is the band dispersion, U the interaction 
strength, a 3 the third Pauli matrix, while the order pa- 
rameter S is defined by (l/N)( £ kQ/3 c l+Qa a lp c k0) with 
N the number of sites. This one-body Hamiltonian can 
be diagonalized by the Bogoliubov transformation 

7ka = "kCka + V k ^(cr 3 ) Q /3C k +Q/3 , 

" (3) 

7ka = U kCka - Uk^(CT 3 ) Q ^C k+Q0 . 



To avoid double counting, k is restricted to the magnetic 
zone. The superscripts c and v refer to the conduction 
and the valence bands split by the exchange Bragg scat- 
tering from the SDW. For simplicity, we assume perfect 
nesting e k = — e k +Q. In this case, the coefficients of 
the transformation are given by it k = [|(f + e k /.E k )], 

v l = [W-^/E*)] > E l = (4 + A2 )> wh ere A = -US/2 
is the SDW gap parameter [9|j. The diagonalized Hamil- 
tonian is given by H = £* Q E * frhatfct ~ TkalL) where 
E k means that the sum extends over the magnetic zone. 
The single-particle energy spectrum is given by ±-E k and 
the SDW ground state for a half-filled band is defined by 
TkcJ^) = TkcJ^) = 0, which corresponds to: 

|fi> = H^KcL - "k^cl^) 1°) ( 4 ) 



= n LK - U k ^2 4+Q/J CT /L c ka)4a 1°) ■ ( 5 ) 



The last form makes the analogies with the BCS ground 
state clear. For example, there exists Andreev-like reflec- 
tions at AF-N interfaces [IcJ • The ground state contains 
coherent particle-hole pairs. To clarify this, we perform 
a particle-hole transformation for states that are in the 
first magnetic-Brillouin zone. Recall that destroying an 
electron in a state, creates a hole in the corresponding 
time-reversed state. For a spinor this can be achieved by 
Ckt ~* ^-kj. an d c ki — > — ^-kf The ground state then 
takes the form 

\ty = n k (w k - u k c k+QT ^L ki )K - uk(i+Qi h -w) \°)h 

(6) 

where c k+ Q Q |0) ft = and /i kQ |0) ft = 0. The particle- 
hole pair is in a triplet state with vanishing net spin pro- 
jection along the quantization axis, has no charge and has 



a wave vector equal to the antiferromagnetic wave vector. 
In the case of a FM, that wave vector would vanish. 

Since we choose a different quantization axis in the two 
AF, we need to include a unitary transformation in spin 
space denoted by U(0, </>) to account for the fact that a t 
spin on one side of the junction is not the same as a f spin 
on the other side of the junction. Since the spin quanti- 
zation axes are taken along the direction of the staggered 
magnetic moment S on each side, the angles (#, <p) corre- 
spond to the orientation of Sr of the right AF expressed 
in the coordinate system of the left side of the junction; 
in Cartesian coordinate S^ = |S^ | (0, 0, 1) and S/j = 
| Sfl| (sin 9 cos (f>, sin 9 sin <j>, cos 9) . The tunneling Hamilto- 
nian then reads H T = {1/N) J] kqro - (lWt kq c k<T d q( 7' + 
/i.e.). The spin flip terms come purely from the choice of 
different quantization axes on the left and on the right. 
There is no real spin flip in the tunneling process. Note 
also that the generators of the unitary transformation 
U(#, 4>) are non-abelian gauge fields that are the ana- 
log of the electromagnetic gauge potentials in the ordi- 
nary Josephson effect. The former do not couple to the 
electromagnetic fields and they do not have a dynamics 
independent of that of the spins. 

Derivation of the equations of motion. The staggered 
magnetic moment operator in the left lead is Si = 
a0 c U+Qa a o'0 Ck where a is a vector of Pauli 
matrices. In the broken symmetry state, the time evo- 
lution of Sl due to Hl is negligible. Since we also have 
[S L ,H R ] = we find dS L /dt = (l/ifi)[S L , H T ], and thus 

dS L /dt = -(V2A0EkqE a/M J^/3 U /55 Wi+QcA<5 - 

/i.e.) whose average S/,(t) = (dS^/di) is given by 

SlO) = ^535Z Im [ <7a / 3U ^* kc l( C k+Qa(*) rf q5(*)) ' 
kq a/35 

(7) 

where (...) is the thermal statistical average with the full 
density matrix. Performing first order perturbation the- 
ory using Ht as the perturbation, one obtains 

/t 
d* , ([4+Qa(*Ki(*).^r(*')]> . 

( 8 ) 

where the average (...} is computed with H = H L + H R 
(the unperturbed part of H). The operators on the right 
are in the interaction representation. 

Once the commutator is evaluated in Eq. J5J, 
one of the factorizations of the four-point correla- 
tion function involves products of correlation func- 
tions on the left and on the right leads such as 
( c k+Qa(*) c ka( t ')}o(rfq+Q5(i)rf q(5 (i')}o- Such correlation 
functions would vanish in a normal paramagnetic state. 
They are non-zero because of the broken symmetry. 
They represent interference in the tunneling process 
between momentum k + Q spin up particles and mo- 
mentum — k spin-down holes, in other words tunnel- 
ing of charge zero spin one S z = coherent particle- 
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hole pairs that have finite momentum and are present 
in the ground-state Eq.Q. In the case of the ordi- 
nary Josephson effect, one would find terms such as 
( c L(*) c -k-<r(*'))o(dq<7(*)d_q_ CT (i'))o that represent tun- 
neling of coherent Cooper pairs. 

In order to compute the averages (...}o in the bro- 
ken symmetry states, we invert the Bogoliubov trans- 
formation Eq. ([3]). Assuming i kq = t kq+ Q = ^k+Q q = 
<k+Qq+Q, we find 

E^q4 +Qc A* = £ E Wr^-M*. (9) 

kq kq ij£{c,v} 



where we defined 

( r kq)cc = ("kWq + <T 3 aa V k Uq + (J^U^Vq + (J 3 aa (J 3 SS V k Vq) 
( r ki)cv = («kUq - CT^ Q W k U q + a%V k Vq - Cr^oJ^ktfq) 
( r kq)l-c = ( W k«q + 0"L U kUq - CT| 5 U k U q - O^CT-fjUkUq) 
( r kq)™ = («k«q ~ CrjL-UkVq ~ ^gVkUq + ^ a O{J«k«q) ■ 

( 10 ) 

Similarly, let fir denote the part of Hr that does not 
commute with c k+ Q Q c? q 5. It can be written as Ht — 

Q-/X) Ekq^' Eij U:^ kq (<) y -7 q WL- Substituting 
these expressions into Eq. ([8|), one finds 



(11) 



where are the Keldysh 



Green func- 
tions in the left (right) lead. Their definitions 

are Sk(q), Q/3 (M') ii^it'hL^M and 

^),c l p^ t ') = - i (7 k (q)a( < )7 k t (q)/5 (^)) I respectively. 
Explicitly, 



analytical result 



G£A* "*) = - /(**)) exp[-i4(t' - *)A]<W . 
- 1) = i/(J& exp[-i^(t' - t)/&]^ , 

where / is the Fermi function. Eq. (jTTJ) for the staggered 
magnetic moment current through a tunnel junction is 
general. A bias could be included. We assume that there 
is no bias so there is no incoherent single-particle tun- 
neling across the antiferromagnetic gap. Integrating over 
t — t' and performing the spin sum in pip , one finds 



where s L(i? ) = S L 



8A L A R 
N 2 



kq 



and 



,/(£ k )-/(-£ k ) 



El) 



(12) 



(13) 



with P indicating principal part. A similar expression is 
found for the equilibrium spin current in the case of fer- 
romagnetic tunnel junctions. Note that the sine function 
present in the standard Josephson case is replaced here 
by a cross product, which is a direct consequence of the 
vectorial nature of the order parameter. 

For a symmetrical junction (A^ = Ar), the same as- 
sumptions and procedure as Ref. [7| lead to the following 



r ir x A(T) tanh(i/?A(T)) 



(14) 



where R = h/ (Aire 2 D 2 \t\ 2 ^j is the (zero-temperature) 

normal-state resistance of the junction with D the den- 
sity of state, which is assumed to be a constant. This 
expression for the temperature dependance of the criti- 
cal current has the same form as that obtained by Ambe- 
gaokar and Baratoff [3] for a BCS superconductor, which 
is not surprising given the formal analogies [ill ]. 

By symmetry, the time derivative of the stagerred mag- 
netic moment on the right lead can be obtained by inter- 
changing the L and R indices in Eq. (p~2|) . As a conse- 
quence, the staggerred magnetic moments of the two AF 
precess about their (constant) sum S^-fS^ at a frequency 
lo = I c \Sl + Sr\/\Sl\\Sr\. In the limit where and 
Sr are nearly colinear, the effect should therefore lead to 
a uniform shift of order u>q in the antiferromagnetic reso- 
nance frequencies of each AF 0. With |^| = U = 2eV 
and a normal state conductance of the order of the con- 
ductance quantum, = 2e 2 //i, one finds a value of 
ujq of order 10 14 Hz, in other words in the visible. This 
frequency would be higher than the single particle gap 
and so would lead to much damping. Resistances that 
are many orders of magnitude larger are thus needed to 
bring its value down. 

AC spin- Josephson effect. In the ordinary Josephson 
effect, the electromagnetic gauge potentials enter directly 
in the argument of the sine function. The present case is 
different. Each magnetic moment associated with a spin 
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couples to the magnetic field through the Zeeman term 
(Hz = — ,9/isB • S) where g is the gyromagnetic ratio 
and \xb the Bohr magneton (in insulators we can neglect 
terms coming from orbital motion). 

Considering magnetic fields Bi and B/? applied re- 
spectively to the left- and right-hand sides of the junc- 
tion, the Heisenberg equations of motion lead to the fol- 
lowing equations of motion for the order parameters : 

Si = -g^B^L x S L + I c s B , x Si, 

Sr = -gfJLB&R x Sfl + I c s L x s R . (15) 

The first term on the right side of the equality is merely 
the contribution of Hz to the Heisenberg equation of 
motion. The second term is the tunneling contribution 
and has exactly the same form as the one we have already 
computed in the zero-field case, namely Eq. (fT2|) and (fl~3|) . 
To see this, return to Eq.©. The only terms from i?r(i') 
that give a non-zero contribution to the average are of 
the form d^ 5 ,(t')ul, a ,(t')c\!_ a '{t'). Since these operators 
are in the interaction representation and [Hq, Hz] = 0, it 
is possible to write d' s ,(t') = <fi s (T')Afg, and Ck<r'(*') = 

A^Ck^r') where = exp[igfj, B B R ^ L ya(t' -^/h]^^ 

and t' stands for the time evolution due to Hq. Sim- 
ilarly, one can show that the unitary transformation 
Ul a , {t') can be written as Kf, 5 ul a {t)h^l, . We then ob- 
tain 4 5 ,(i04 CT '(*')cw(*') = dl s (r')ul(t)c^(r'), U 
all spin indices are at time t. Therefore, it is possible to 
factor the unitary transformation U $ (t) out of the inte- 
gral over t', so that the rest of the calculation just follows 
the same path as in the zero-field case, thus leading to 
Eq. (HU). 

If Bi = Bfl = B, Eq. (fr5|) implies that in the rotating 
frame u = {x',y',z'} defined by du/dt = — g/i^B x u, 
Sl and Sr still precess about their (constant) sum S = 
Sl + S/j at a frequency ujq. Returning to the static frame, 
Eq. (fT5|) gives dT,/dt = -gfi B & x X so that Si and S^ 
undergo a motion of double precession. 

In the ordinary Josephson effect, a constant electric 
potential difference V leads to a phase difference ip that 
depends linearly on time (ip — 2eV/K). In the present 
case, an analog is found by computing the time depen- 
dence of the relative orientation $ between Si and Sr. 
Using Eq. (fTS"]) to compute d(cos-d)/dt = e?(§i • s R )/dt, 
we obtain 

= -g/i B SB ■ e(i), (16) 

where <5B = Bi — B/? and e(t) = Si(t) x Sfl(t). As in 
the ordinary Josephson effect, this expression does not 
contain cxplicitely the tunneling matrix element. Con- 
trary to the ordinary Josephson effect, this equation is 
non-linear. Solving it numerically along with Eq. (|15p 
in the case where the magnetic field vanishes on one side 
of the junction, one finds that the angle between Si and 
Sr behaves as a sine-like function of time. The presence 



of an additionnal constant magnetic field B throughout 
the system adds a beat to this sine-like behavior. The 
gyromagnetic ratio g is material dependent. No such 
non-universal constant appears in the ordinary Joseph- 
son effect. The above discussion of the AC effect can be 
transposed for FM by replacing the staggered magnetic 
moment by the uniform one. 

In summary, we found a Josephson-like equilibrium 
current of the staggered magnetic moment through a tun- 
nel junction between two itinerant AF with non-colinear 
staggered magnetic moments. In analogy with the fer- 
romagnetic case, this current is proportional to Sr x Si 
where Si and Sr are the staggered magnetic moments 
on either sides of the junction. The effect comes from 
tunneling of coherent spin-one neutral particle-hole pairs 
with a net momentum and no net spin projection in the 
order parameter direction. There are important differ- 
ences with the ordinary Josephson effect coming both 
from the different coupling of the electromagnetic fields 
and from the non-abelian nature of the broken symmetry. 
It would be interesting to perform magnetic resonnancc 
experiments to detect the effects we predicted. 
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